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We formulate a theory of non-Abelian fractional quantum Hall states by considering an anisotropic 
system consisting of coupled, interacting one dimensional wires. We show that Abelian bosoniza- 
tion provides a simple framework for characterizing the Moore Read state, as well as the more 
general Read Rezayi sequence. This coupled wire construction provides a solvable Hamiltonian 
formulated in terms of electronic degrees of freedom, and provides a direct route to characterizing 
the quasiparticles and edge states in terms of conformal field theory. This construction leads to a 
simple interpretation of the coset construction of conformal field theory, which is a powerful method 
for describing non Abelian states. In the present context, the coset construction arises when the 
original chiral modes are fractionalized into coset sectors, and the different sectors acquire energy 
gaps due to coupling in "different directions" . The coupled wire construction can also can be used 
to describe anisotropic lattice systems, and provides a starting point for models of fractional and 
non-Abelian Chern insulators. This paper also includes an extended introduction to the coupled 
wire construction for Abelian quantum Hall states, which was introduced earlier. 
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I. INTRODUCTION 

The search for non-Abelian states in electronic mate- 
rials is an exciting frontier in condensed matter physics 1 . 
Motivation for this search is provided by Kitaev's 
proposaP to use such states for topological quantum 
computation. The quantum Hall effect is a promising 
venue for non Abelian states. There is growing evi- 
dence that the pfafHan state introduced by Moore and 
ReacpE' describes the quantum Hall plateau observed 
at filling v = 5/#®. The Moore Read state gives the 
simplest non-Abelian state, with quasiparticles that ex- 
hibit Ising non-Abelian statistics. While the observa- 
tion and manipulation of Ising anyons is an important 
goal, Ising anyons are not sufficient for universal quan- 
tum computation^. The Z% parafermion state introduced 
by Read and Rezayf^ is a candidate for the quantum 
Hall plateau at v = 12/5. The quasiparticles of th e Read 
Rezayi state are related to Fibonacci anyons^El which 
have a more intricate structur e tha t in principle allows 
universal quantum computation^^. 

There is currently great interest in realizing quantum 
Hall physics in materials without an external magnetic 
field or Landau levels. This possibility was inspired by 
Haldane's realizatiorfS' that a zero field integer quantum 
Hall effect can occur in graphcnc, provided time reversal 
symmetry is broken. Though such an anomalous quan- 
tum Hall effect has not yet been obser ved, related physics 
occurs in topological insulator d 15 * 16 ! , which have been 
predicted and observed in both two and three dimen- 
sional systems. Recently, there have been suggestions 
for generalizations of this idea to zero field fractional 
quantum Hall states^nUl as well as fractional topological 
insulators^Jl^iJ. The question naturally arises whether it 
is possible to e ngineer zero field non-Abelian quantum 
Hall statesP 1221 , which exhibit the full Ising, or even Fi- 
bonacci anyons. 

A difficulty with answering this question is a lack of 



methods for dealing with strongly interacting systems. 
Moore and Read, and later Read and Rezayi, built on 
Laughlin's idea^ and constructed trial many body wave 
functio ns a s correlators of non trivial conformal field 
theories^EHl. This allowed most properties of the quasi- 
particles and edge states to be deduced, and it had the 
virtue of allowing the construction of interacting elec- 
tron Hamiltonians with the desired ground state. How- 
ever, since this approach relied on the structure of the 
lowest Landau level, it is not clear how it can be ap- 
plied to a lattice system. Effective topolo gica l quantum 
field theories^HSl anc i parton constructions^ provide an 
elegant framework for classifying quantum Hall states 
and provide a description of their low energy properties. 
However, since the original electronic degrees of freedom 
are replaced by more abstract variables, these theories 
provide little guidance for what kind of electronic Hamil- 
tonian can lead to a given state. 

In this paper we introduce a new method for describ- 
ing non Abelian quantum Hall states by considering an 
anisotropic system consisting of an array of coupled one 
dimensional wires. Study of the anisotropic limit of quan- 
tum Hall states dates back to Thouless et alP2 w ho used 
this limit to evaluate the Chern invariant in the inte ger 
quantum Hall effect. The Chalkcr Coddington model 35 
for the integer quantum Hall effect also has a simple 
anisotropic limit, which is closely related to the coupled 
wire model. A coupled wire construction for Abelian frac- 
tional quantum Hall states was introduced in Ref. 1361 
Here we build on that work and show how the coupled 
wire construction can be adapted to describe the Moore 
Read pfaffian state, as well as the more general Read 
Rezayi sequence of quantum Hall states. 

The coupled wire construction has a number of de- 
sirable features. First, it allows for the definition of a 
simple Hamiltonian, expressed in terms of electronic de- 
grees of freedom, that can be transformed, via Abelian 
bosonization, into a form that for certain special param- 
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eters can be solved essentially exactly. The quasiparticle 
spectrum, as well as the edge state structure follow in 
a straightforward manner. The coupled wire construc- 
tion thus provides a direct link between a microscopic 
electron Hamiltonian and the low energy conformal field 
theory description of the edge states. As such, it provides 
an intermediate between the wave function approach to 
quantum Hall states and the effective field theory ap- 
proach. 

The coupled wire construction also provides a simple 
picture for the quantum entanglement present in quan- 
tum Hall states^HSS When the electron Hamiltonian is 
transformed via Abelian bosonization, it becomes iden- 
tical to a theory of strips of quantum Hall fluid cou- 
pled via electron tunneling between their edge states. It 
thus provides a concrete setting for the more abstract 
coupled edge state models considered by Gils, et alP^. 
The coupled wire model is also similar in spirit to the 
AKLT model of quantum spin chains^, which provides 
a similarly intuitive and solvable model for understanding 
fractionalization in one dimension. For the non Abelian 
quantum Hall states, the coupled wire construction pro- 
vides a concrete interpretation for the coset construction, 
which is a powerful (albeit abstract) mathematical tool 
for describing non Abelian quantum Hall states^!. 

A final virtue of the coupled wire construction is that 
it can be applied to zero field anisotropic lattice models. 
The effect of the magnetic field in the coupled wire model 
is to modify the momentum conservation relations when 
electrons tunnel between wires. A similar effect could 
arise due to scattering from a periodic potential. The 
coupled wire construction may thus provide some guid- 
ance for the construction of lattice models for the frac- 
tional and non Abelian quantum Hall states. Note that 
our construction is somewhat different from the propos- 
als for Chern insulators Refs. 17-20| because we do not 
require nearly flat bands with an non zero Chern number. 

The outline of the paper is as follows. We will begin 
in section II with an extended introduction to the cou- 
pled wire construction for Abelian quantum Hall states. 
Much of this material was contained either explicitly or 
implicitly in Ref. 1361 Here, since we are free from the 
constraints of a short paper, we will fill in some details 
that were absent in Ref. 1361 In particular, we will ex- 
plain the generalization of the coupled wire construction 
to describe systems of bosons, we will demonstrate the 
Abelian fractional statistics of quasiparticles described 
in our approach, and we will explicitly construct second 
level hierarchical fractional quantum Hall states. 

Section III is devoted to the Moore Read state. We 
will begin with a construction of this state for bosons 
at filling v = 1. This leads to a bosonized model that 
can be solved using via fermionization. We will then 
show that for a special set of parameters the model has a 
particularly simple form, which can be interpreted in the 
framework of the coset construction of conformal field 
theory. We will conclude section III by showing how to 
construct the more general Moore Read state at filling 



v = 1/(1 + q) for q even (odd) for bosons (fermions). 

In section IV we will generalize our construction to 
describe the Read Rezayi sequence at level k. Again, 
the formulation is simplest for bosons at filling v = k/2, 
where the coupled wire model is closely related to the 
coset construction of these state. We show that the cou- 
pled wire model leads maps to a bosonized representation 
of the critical point of a Zk statistical mechanics model, 
which for fc = 3 reduces to the 3 state Potts model. 
This bosonized representation allows us to identify the 
Zk parafermion primary fields and fully characterize the 
edge states of the Read Rezayi states. Finally, as in sec- 
tion II, we conclude by generalizing our results to de- 
scribe bosonic (fermionic) level k Read Rezayi states at 
filling k/(2 + qk) for q even (odd). 

Some of the technical details are presented in the ap- 
pendices. Appendix A gives a careful treatment of Klein 
factors, while Appendices B and C contain some of the 
conformal field theory calculations discussed in section 
IV. 



II. ABELIAN QUANTUM HALL STATES 

A. Coupled Wire Construction for Fermions 

In this section we review the coupled wire construc- 
tion for fermions introduced in Ref. 36, We begin by 
considering an array of identical uncoupled spinless non 
interacting one dimensional wires, as shown in Fig. la, 
with a single particle electronic dispersion E(k). We as- 
sume each wire is filled to the same density, characterized 
by Fermi momentum kp. The two dimensional electron 
density is then n e = kp/i^a, where a is the separation be- 
tween wires. A perpendicular magnetic field, represented 
in the Landau gauge A = ~Byx, shifts the momentum 
of each wire. The right and left moving Fermi momenta 
of the j'th wire are then 

k*( L =±k° F + bj, (2.1) 

where b = \e\aB/h. The filling factor hn e /\e\B is then 

v = 2k F /b. (2.2) 

The low energy Hamiltonian, linearized about the Fermi 
momenta is 

Ho = j dxv °F (iP?H-id x - fc* )Vf 

-i>f L \-id x -k^)rf), (2.3) 

where ipj L describe the fermion modes of the j'th wire 

in the vicinity of the Fermi points kpj L . 

We next bosonize by introducing bosonic fields <fj{x) 
and 0j(x) that satisfy 

[d x 6j(x),cpj'(x')] =i%5jf8 xx >, (2.4) 



where we use the shorthand notation 5 XX ' = 8(x — x'). 
(pj(x) is a bosonic phase field, while 9(x) describes den- 
sity fluctuations. The long wavelength density fluctua- 
tions on the j'th wire are 

Pj (x) =Y,i>f{xWM = dJiix)/*. (2.5) 

a 

The electron creation and annihilation operators may be 
written 



y/2nx c 

3 y/2lTX c 



(2.6) 



where x c is a regularization dependent short distance cut- 
off and Kin IS El Klein factor that assures the anticommuta- 



tion of the fermion operators on different wires. Eq. 2.4 
hides the zero momentum parts of 9j and <pj , which must 
be accounted for in order to correctly treat the Klein fac- 
tors. Since this issue tends to obscure the simplicity of 
our construction, we will not dwell on it in the text of 
the paper. Appendix A contains a careful discussion of 
the zero modes and Klein factors, which shows when they 
can be safely ignored. 

In terms of the density and phase variables, the Hamil- 
tonian for non interacting electrons is 



V. 



2tt ^ 

3 



dx [(d x ipj) 2 + (d x t 



(2.7) 



Interactions between electrons as well as electron tun- 
neling between the wires can be added. In general, there 
are two classes of terms : forward scattering and inter 
channel scattering. The forward scattering terms con- 
serve the number of electrons in each channel and can be 
expressed as the interactions between densities and cur- 
rents. This leads to a Hamiltonian that is quadratic in 
the boson variables, 



E 

jk 



J dx ( d x (pj <),(), ) M jk 



(2.8) 

Here the 2x2 matrix 'M.jk = Sjkivp/2Tr + Ujk, where 
TJjk describes the forward scattering interactions. Hg LL 
describes a gapless anisotropic conductor in a sliding Lut- 
tinger liquid phased®'. 

Symmetry allowed inter channel scattering terms must 
be added to "Hg LL . They can open a gap and lead to inter- 
esting phases. The allowed terms are built from products 
of single electron operators, and have the form 



(2.9) 



where Sp^ L are integers such that ipf+p (ipf+£^) appears 
|s^ L | times for Sp^ L > (< 0). It is convenient to write 



(a) 



(b) 








K F,j K F,j 



K F.j+1 K F.j+1 



FIG. 1. (a) An array of coupled wires in a perpendicular mag- 
netic field, (b) The magnetic field shifts the momentum of the 
wires. At at special filling factors, there exist momentum con- 
serving correlated tunneling processes that lead to quantum 
Hall states. The process shown describes the Laughlin state 
at filling v = 1/3. 



Sp'^ L in terms of a new set of integers, 
Sp = (rip + mp)/ 2 
s p = ( u p ~ m p)/ 2 - 



Then o\ mp ' np} {x) takes the form 



{m P ,n p } = cRjRexpi 



(2.10) 
(2.11) 



(2.12) 



where c is a non universal constant. The product of Klein 
factors is 



~{m p ,n p } = -J-j- ^ 



j+V 



(2.13) 



The oscillatory factor describing the net momentum of 



,{m p ,n p } . 



R{ m »< n *y ( x ) = expi b P n P + k° F m p ^j 



x. (2.14) 



The operators £)j m^ " n3, } define a term in the Hamiltonian, 



o-n P } = J2j dx (v {m ^ } e>j mp '" p} (x) + h. 



(2.15) 

There are physical constraints on the allowed {m p , n p }. 
Since Sp^ L must be integers, we require 

Trip = rip mod 2. (2-16) 

Charge conservation requires that 

^% p = 0. (2.17) 
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Momentum conservation implies 

y^(bpn p + k F m p ) = 0, 



(2.18) 



so that the oscillatory term in (2.14) vanishes. 
The Hamiltonian 



{m p ,n p } 



(2.19) 



can be studied perturbatively using the standard renor- 
malization group analysis. The lowest order RG flow 
equation for v^ mp ' np ^ is 



^{TOp.np} /(j£ = (2 — A^ mp ' np ^ )v^ mp ' np ^ . 



(2.20) 



The scaling dimension A^ mp/np ^ depends on the forward 
scattering interactions M.jk in Hg LL . When A^ mp ' np ^ > 
2 the operator 0^ m " np ^ is irrelevant and v^ mp - n "^ does 
not destabilize the gapless sliding Luttinger liquid fixed 
point. When A^ mp ' rlp ^ < 2, 0l m p>™p} is relevant and 
v {m p ,n p } g rows a ^ \ ow energy, destabilizing the sliding 
Luttinger liquid. In principle rVL-j. can be parameter- 
ized given an underlying model of the electron-electron 
interactions. However, 'M.jk may be renormalized by ir- 
relevant and/or momentum non conserving operators, so 
it may not resemble the bare interactions. Here we fol- 
low the approach of Ref. [35] and assume that Mjfe have 
values such that a particular operator (or set of opera- 
tors) 0{™p'™p} is relevant. Our object is to characterize 
the resulting non trivial strong coupling phases. There 
are special values of Mj/- that lead to particularly simple 
boson Hamiltonians that can be solved exactly. These 
solvable points provide a powerful way to characterize 
the resulting strong coupling phases. 

As shown in Ref. (35] a number of non trivial 2D 
phases can be analyzed using this approach, including 
Abelian fractional quantum Hall states, superconductors 
and crystals of electrons, quasiparticles or vortices. In 
particular, Abelian quantum Hall states are described by 
a single relevant operator {m p , n p } satisfying J2 P n p ^ 0. 
From (2.2) and (2.18) this corresponds to a filling factor 



v = 2 



(2.21) 



In Section II. C we will review this construction for the 
Laughlin states and the Abelian hierarchy states. But 
first, we will show that the coupled wire construction can 
also be straightforwardly applied to systems of bosons. 



density variable 0{x) satisfying (2.4). The density flue 



tuations have important contributions near wavevectors 
q n = 2imp that are multiples of the average ID density 



p(x) = p + ^2p n (a 



(2.22) 



As with the fermions, the long wavelength density fluctu- 
ation is po(x) = d x 9(x)/ir. The density wave at q ~ 2irpn 
is, 



(2.23) 



Here and in the following we will denote the ID density 
p in terms of "2fcp" = 27rp. This allows us to proceed 
analogously with the fermions and use formulas (2.2 ) and 
pT8| for the filling factor. 



The Hamiltonian for bosonic wires coupled only by 
long wavelength interactions has exactly the same form 
as "Hg LL . The only difference is that the non interacting 



Pauli compressibility term Hz. 



is absent. Tunneling 



a boson between wire j and j + p in the presence of a 
magnetic field is described by the operator 



^ j+p (x)^(xy b ^, 



(2.24) 



where <&j(x) oc expiipj(x) is the boson creation operator. 
Due to interactions this process can involve scattering 
from the 2kpn density fluctuations of the bosons. The 
most general coupling term thus has the form 



O 



= cR {r 



exp i 



IpOj+p I ) 



where R is given in (2.14) 



(2.25) 

This is almost identical to 
the inter channel scattering terms for fermions. The only 
differences are the absence of Klein factors and the con- 
straints on the allowed values of {m p , n p }. Charge and 
momentum conservation still requires (2.17 2.18), but 
unlike for fermions, where m p and n p 
corresponding constraint for bosons is 



obey (2.16), the 



= mod 2. 



(2.26) 



The analysis of bosonic states then follows in exactly the 
same manner as fermionic states, as described in Eqs. 



(2.19)-(2.21). 



C. Laughlin States v = 1/m 



B. Coupled Wire Construction for Bosons 

We now consider coupled wires of one dimensional 
bosons. The low energy excitations of a single wire can 
be described by "bosonizing the bosons" , to express them 
in terms of a slowly varying phase (p(x) and a conjugate 



Here we will examine the coupled wire construction for 
the Laughlin states in some detail. We include the details 
here because the Laughlin states provide the simplest non 
trivial application of the coupled wire construction. We 
begin by introducing the relevant interaction term, and 
then characterize the bulk quasiparticles and edge states. 
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(a) 





Oj+1/2 1 

^ 



(b) 
j + 1 



o 



j-1/2 



j-1 



Oj+1/2 


v=1/m 


: i : 








v=1/m 


z i r 


Oj-1/2 


v=1/m 



• f-1 



FIG. 2. (a) Schematic representation of the correlated tunnel- 
ing processes in (2.27 1 that lead to the Laughlin state. Verti- 
cal arrows represent <f> a nd cir cular arrows represent 9. After 
the change of variables (2.28 1, the model describes strips of 
v = 1/m quantum Hall fluid coupled by tunneling electrons 
between their edge states. 



1. Tunneling Hamiltonian 

The Laughlin sequence of quantum Hall states at filling 
v = 1/m are characterized by the correlated tunneling 
operators involving two neighboring wires. The relevant 
operator is associated with a link t = j + 1/2 between 
wires j and j+1, 



Oe = expi [ipj — ipj+i + m(9j + 9 J+1 )] 



(2.27) 



Using (2.21) it can readily be seen that this term is 
allowed for magnetic fields b corresponding to v = 1/m. 
Moreover, from (2.161 and (2.26), it is clear that m odd 
(even) corresponds to a fermionic (bosonic) state, as is 
expected for the Laughlin state. In Eq. |2.27| we have 
suppressed the Klein factors, which are necessary for 
fermions. They are treated carefully in Appendix A. This 
term is represented schematically in Fig. 2a. The no- 
tation for this diagram is slightly different from the one 
used in Ref 36, The vertical arrows describe the tunneling 
of charge between the wires (represented by <frj+i — <frj), 
while the circular arrows describe backscattering within 
a wire (represented b y fl ? -j + i). Note that the number of 



#'s is constrained by (2.16) or (2.26). 

The tunneling operators defined above have the spe- 
cial property that they all commute with eachother. In 



particular, [O 



2] = 0. This means that the 



i + l/2,Cj-l/2 

components from wire j in those two operators commute 
with one another. This invites us to introduce right and 
left moving chiral fields on wire j that distinguish the 
two contributions to 0j±i/2- We thus write 



bj = <pj + m6j 



ipj — mOj 



if 



(2.28) 



The decoupling can be explicitly seen from the commu- 
tation algebra, 



d x 4%(x),<t%, (x r ) = 2-Kimp8 pp >5jj>8 xx <- 



(2.29) 



The interaction term is now 0j + i/2 = expi(^ +1 — cj)^). 
The charge density is pj = {d x ^)f — d x M) / {2-Km) . 

It is also convenient to introduce new density and 
phase variables defined on the links I = j + 1/2 between 
wires, 



§t = (4>f - # +1 )/2 



(2.30) 



These satisfy [d x 9 e (x),9 e (x')} = [d x {p t (x),(p v (x')] = 
and 



[d x 9i(x),tpt'(x')] = inm5 U i5 x 



(2.31) 



The charge density associated with the link I can be writ- 
ten 



p e = d x 9~i/(mir). 



(2.32) 



In terms of the new variables, the Hamiltonian be- 
comes, 

H = Hsu J [9,0\+^2 f dxvcos29 e , (2.33) 

t 

where without loss of generality we have assumed v is 
real. As shown in Appendix A, there is no Klein factor, 
provided the zero momentum component of 9g is correctly 
defined. 

Provided the forward scattering interactions defining 
^SLL are sucn that v is relevant, the system will flow at 
low energy to a gapped phase in which 9g is localized in 
a well of the cosine potential. As argued in Ref. [351 it 
is always possible to find such interactions. In particu- 
lar, consider a simple interaction such that Hsll has the 
decoupled form 



Hi 



SLL 



«0 \ - 
2tt 2-« 



dx 



-(dJif 
9 



g(d x tpe 



(2.34) 



The scaling dimension of cos 29g is A = mg. It follows 
that for g < 2/m, v is relevant. It should be empha- 
sized that "Hg LL can be expressed in terms of the original 
fermion operators, which includes a specific four fermion 
forward scattering interaction. For special values of g 
this model can be solved exactly. In the limit g — > 0, the 
variable 9i becomes a stiff classical variable, so that the 
approximation of replacing — cos 29 1 by 29f becomes ex- 
act. For larger g, we rely on our understanding that g is 
renormalized downward by v, so that at 9 stiffens at low 
energy. For g = 1/m there is another exact solution be- 
cause it is possible to define new variables such that the 
Hamiltonian has precisely the form of (2.7). The prob- 



lem can then be refermionized and expressed in terms 
of non interacting fermions which have a single particle 
energy gap. We will not dwell on these exact solutions 
any further in this paper. We will be content with our 
understanding that any g < 2/m leads to a g apped state. 

The gapped phase is the Laughlin stated This can be 
seen by examining the quasiparticle excitations and the 
edge states. 



6 



2. Bulk Quasiparticles 

Quasiparticles occur when Qi{x) has a kink where it 
jumps by n. From (2.32) it can be seen that such a kink 
is associated with a charge e/m. This makes the charge 
fractionalization in the fractional quantum Hall effect ap- 
pear similar to the fractionalization that occurs in the one 
dimensional Su Schrieffer Heeger (SSH) model 49 . How- 
ever, there is a fundamental difference. The solitons in 
the SSH model occur at domain walls separating physi- 
cally distinct states. This prevents solitons from hopping 
between wires via a local operator. In contrast, the states 
characterized by 9g and 9i + 7r are physically equivalent. 
They are related by a gauge transformation in which, say, 
ifj — > (fj + 27r, which takes to 9j±i/2 T 71 '■ This 

allows quasiparticles to hop via a local operator without 
the nonlocal string. Though SSH solitons and Laugh- 
lin quasiparticles are distinct, they become equivalent on 
a cylinder with finite circumference. The Tao Thouless 
"thin torus" limitP^ can be described the extreme case in 
which the "cylinder" consists of a single wire with elec- 
tron tunneling "around" the cylinder. In this case, our 
theory maps precisely to an m state version of the SSH 
model. 

The local operator that hops quasiparticles between 
links j + 1/2 and j — 1/2 is simply the backscattering of a 
bare electron on wire j, ipjip^ (or equivalently for bosons 
the 2kp density operator). Using the transformations 



To move a quasiparticle from x\ to xi on link t\ and 
then to Xi on £2, use the operator 



£2-1/2 

Pl 1 {xi,X 2 ) Yi Xj{x2)- 

Since 9 is gapped, this can be written 



(2.38) 



2i(e e (x 2 ))/m 



(2.39) 



The string of (9) is responsible for the fractional statis- 
tics. 

Consider moving a quasiparticle through a closed loop. 
The operator that takes a quasiparticle around the rect- 
angle for med b y Xi, x%, l\ and £2 can be constructed 
by doing (2.38) twice, which eliminates the quasiparticle 



operators. This then gives a phase 
< 2 -i 



n 

i=ii+i 



e 2i({e e (x 2 ))-{e e {x 1 )))/m _ e 2niN QP /m 



(2.40) 



(2.28) and (2.30) it is straightforward to show that 



Xj [x') = tf^f 



e 2i8 3 (x') = e i(4>f(x')-^(x'))/ m 
e i(<Pj + l/2-<Pj-l/2+8j + l/2+8 1 -l/2)/m_ 



(2.35) 



where Nqp is the number of quasiparticles enclosed by 
the rectangle. Here we have used the fact that ((9g(xx)} — 
(9t(x2)))/2TT simply counts the number of quasiparticles 
on link £ between x\ and X2- 



3. Edge States 

For a finite array of wires with open boundary condi- 
tions, the edge states are apparent, since there are extra 
chiral modes left over on the first and last wire. From 



From (2.31) it can be seen that this operator takes (2.29), it can be seen that these modes have precisely 



d x 6j±x/2 to d x 9j±i/2T^S(x— x'), transferring a quasipar- 
ticle from j + 1/2 to j — 1/2. The operator that transfers 
a quasiparticle along wire from x to x' along wire link £ 



the chiral Luttinger liquid structure of v — 1/m edge 
statesP. 



H 



edge 



mVQ 
47T 



(dx^i 



L\2 



(2.41) 



Pi{x,x') 



i((pi(x)-(pe(x')) /m _ i J*i dxd x ipt/m 



(2.36) 



which can also be expressed in terms of the bare electron 
densities and currents. 

A quasiparticle operator may be defined as 



with [d x <P?(x),cj>X(x')} 
on the j = 1 edge is 



2iri8 xx i . The electron operator 



(2.42) 



R/L 

QPJ+1/2 



(X) 



1 R / L 1 

p . ... Ti 



i(<fi j+1/2 +/-e j + 1/2 )/m 



(2.37) 



It is straightforward to show that this operator has the 
expected dimension A — m/2, characteristic of the chiral 
Luttinger liquid. 



In the bulk, since 9 is gapped, we have ^qp ( 

r/l 



One can view the change of variables (2.28) as a trans- 



QP,i 



3 2i(0c) _ of coursej s i nC e ^>Q Pe can not be locally 



built out of bare electron operators, it is not by itself a 
physical operator. However, the operator that transfers 
a quasiparticle from one location to another can be built 
from a string of local operators like (2.35) and (2.36). 



This allows the fractional statistics of the quasiparticles 
to be seen quite simply. 



formation between a sliding Luttinger liquid built out of 
bare electrons and a sliding Luttinger liquid built out of 
v = 1/m edge states. The correlated tunneling term for 
the bare electrons becomes the electron tunneling oper- 
ator coupling the edge states. The array of wires then 
becomes an array of strips of v = 1/m quantum Hall 
fluid coupled by electron tunneling, as shown in Fig. [2] 
When the electron tunneling is relevant the strips merge 
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FIG. 3. (a) Schematic of tunneling processes in (2.451 that 
lead to 2nd level Abelian hierarchy states, (b) After the 
transformation ( |2.47[ ) the model describes coupled strips of 
v — In /(mo + mi) quantum Hall state coupled by tunneling 
electrons between the two channels of edge states. 



to form a single bulk v 
modes at the edge. 

The quasiparticle operator at the j 



1 jm fluid, leaving behind chiral 
1 edge is 

(2.43) 



*QP,1 — e 



As discussed above, since ^qp can not be made out of 
bare electron operators, it is not by itself a physical op- 
erator. However, quasiparticle tunneling from the top to 
the bottom edge can be built from a string of backseat - 
tering operators (2.35). When the gapped bulk degrees 
of freedom are integrated out, this string of operators 
becomes 



N 

4=1 



it 



QP,l^QP,N- 



(2.44) 



D. Hierarchy States 

In this section we show how the coupled wire construc- 
tion descr ibes h ierarchical Abelian fractional quantum 
Hall states^ES We restrict ourselves to second level 
states, which are characterized by a 2 x 2 K matrix^. 
Generalization to higher levels is straightforward. 

2nd level hierarchy states arise from an interaction 
term that involve three coupled wires. A generic term 
is shown in Fig. [3j and can be described by the operator 

Oj = expi [n(<pj-i — tpj+i) + 2m Bj + m 1 (9 J+ i + . 

(2.45) 

Here n and mo are any integers, while mi (mi + n) is 
an even integer for bosons (fermions). Again, we defer 
discussion of the Klein factors to Appendix A. From[2~2~l| 



this interaction conserves momentum at a filling factor 

2n 



m 



(2.46) 



This set of states corresponds to the standard Haldane- 
Halperin hierarchy states at filling (po + 1/pi)^ 1 (po is 
even (odd) for bosons (fermions) and p\ is even) for the 
choice, n = px/2, mo = poPi/2 and mi = popi/2 + 1. 

To analyze this state, we group the wires into pairs 
j = 2k and j = 2k + 1. Pair k is connected to pair k + 1 
by two tunneling terms, 02/c and 02/c+i- As in (2.28) we 



define new variables that decouple right and left moving 
modes on the pairs of wires. 



2m 



<j) kl = nip 2 k-i + mi# 2 fc-i 
4>k,i = nip 2 k-i ~ mi6» 2 fc-i 
4>k,2 = n( P2k + ■mi9 2 k 
<f>k,2 = nL P2k - mi6 2 k - 2m 6 2 k-i- 
The new fields obey the commutation algebra 

9x4>l ia (x),(i) p k , b {x') = 2irip5 pp i5 k k>K ab 8 xx , 
where the K matrix is 
K ab = 
The charge density is 

[t a d x (^ M -4, a )/2n 



(2.47) 



mi m 
mo m i 



Pk 



with 



tn 



1 



m + mi 

We next define variables on links i = k + 1/2, 



t>a,e=k+l/2 
¥>a,e=k+l/2 



6£. 



0fc+l,a)/2 



(2.48) 



(2.49) 



(2.50) 



(2.51) 



(2.52) 
(2.53) 



These satisfy, [d x La ,Q t > .b] = [d x <Pl,a,<Pi\b[ = and 

[d x 6i >a ,tpt> tb ] = iirK ab 6w>5 XX ' . (2.54) 

In terms of these new variables, the Hamiltonian in the 
presence of the correlated n-electron tunneling operators 
becomes 



E 



dxv ( cos 29& i + cos 28. 



t:2 



(2.55) 

If v flows to strong coupling, we have a gapped bulk, 
describing a v — 2n/(mo + mi) quantum Hall fluid char- 
acterized by the K matrix ( |2.49 1. As in Section II. C. 3, 
this can be interpreted as quantum Hall strips with edge 
states coupled by the charge ne tunneling operators 



ne,Rf L 
k,a 



R/L 



(2.56) 
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Quasiparticles, given by ir kinks in Oi i or 9i 2 are created 
by 



y QP,oi+l/2 ~~ e 



(2.57) 



They have charge e/(mo + mi). The bare electron 
backscattering operator corresponds to quasiparticle tun- 
neling, 

Xk, a = e 2 ^-^ = *^ a , fc _ 1/2 *§P )0)fc+1/2 . (2.58) 

At the edge of a semi infinite system, where will be two 
chiral modes left over described by (j> R a - From (2.48) it 
can be seen that these give precisely the chiral Luttinger 
liquid edge states characterized by the K matrix (2.49). 



III. MOORE READ STATE 

We now generalize the coupled wire construction to 
describe the Moore Read stated Our approach was 
motivated by the observation by Fradkin, Nayak and 
Shoutens^ that the Moore Read state for bosons at fill- 
ing v — 1 has a simple interpretation in terms of two 
coupled copies of bosons at v = 1/2. Each copy is de- 
scribed by a SU(2)i Chern Simons theory, and the cou- 
pling between them introduces the symmetry breaking 
SU(2) 1 x SU(2) 1 -> SU(2)^. 

We therefore first consider the problem of coupled 
wires of bosons at filling v = 1, where the bosons on 
each wire have two flavors, each at v = 1/2. The allowed 
boson tunneling and backscattering terms in our con- 
struction have a simple representation in the low energy 
bosonized theory. Moreover, by fermionizing the bosons, 
the Majorana fermions associated with the Moore Read 
state^ emerge naturally. 

There is a special set of values for the interactions in 
which the problem is particularly simple. In this case, the 
Hilbert space associated with the two right (left) moving 
chiral modes on each wire decouples into two sectors. 
One of the sectors is coupled to the corresponding sector 
of the left (right) moving modes on the same wire, while 
the other sector is coupled to the corresponding sector 
of the left (right) moving modes on the neighboring wire. 
Both of these couplings introduce gaps, but the two sec- 
tors are gapped in "opposite directions". This gives a 
kind of hybrid between the insulating phase, in which all 
chiral modes are paired on the same wire, and the quan- 
tum Hall states, in which all the chiral modes are coupled 
on neighboring wires. What gets left behind on the edge 
is a fraction of the original chiral modes. 

This fractionalization of the original chiral modes de- 
scribed mathematically i n term s of the coset construction 
in conformal field theory^EH. The original pair of chi- 
ral modes are described by a SU(2)\ x SU{2)i theory 
with central charge c = 2. These modes decompose into 
three sectors: SU(2) 1 x SU(2) 1 / SU(2) 2 , SU(2) 2 /U(l) 
and U(l) with c = 1/2, 1/2 and 1, respectively. The inde- 
pendent sectors are then gapped "in different directions" . 



We will describe this construction in Section III.B. This, 
in effect, gives a concrete and somewhat more explicit im- 
plementation of the Fradkin, Nayak, Shoutens construc- 
tion. 

After establishing the Moore Read state for v = 1 
bosons, we will go on to generalize our construction to 
account for fermions, and the q-pfaffian state at filling 
v = l/(q+l), where q is even (odd) for bosons (fermions). 



A. Bosons at v = 1 

We begin with a Hamiltonian H = Hg^+V describing 
coupled wires of two component bosons, which can be 
viewed as a double layer system, as in Fig. Each 
component has a density corresponding to filling v = 1/2. 

except now each wire has 



H° has the same form as (2 



two bosons, 0j^ a and (fj t a, for a = 1, 2, which satisfy 

[d x j: a(x),(p 3 ',a'(x')} = iTTdjj'Saa'Sxx'- (3.1) 

The interaction terms V consist of boson tunneling and 
backscattering operators that are consistent with mo- 
mentum conservation. We consider three such terms, 
depicted in Fig. [4] 




tabOj t ab 



h.c 



(3.2) 

The first term involves coupling between channel a on 
one wire and channel b on the neighboring wire. 



j;ab 



ai(<Pj,a—<Pj+l,b+2( 



!>)) 



(3.3) 



This term is similar to (2.27). The coefficient 2 of the 9 



terms is fixed by the filling factor. In addition, there are 
allowed terms that couple the two channels on a single 
wire. These include a Josephson like coupling between 
the two channels, 



qu _ e i(tpj,i—(pj,a) ^ 



(3.4) 



as well as an interaction that locks the u 2kp" densities 
of the two channels, 



qv _ e i(%,i-%,2) 



(3.5) 



These three terms (as well as combinations of them) are 
the only allowed interaction terms at v = 1/2 that in- 
clude up to first neighbor coupling. 

It is now useful to introduce right and left moving chi- 
ral fields, 



26 



26>., 



<Pj,a 

as well as "charge" and "spin" fields, 

€ 2 )/2- 



(3.6) 



(3.7) 
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FIG. 4. (a) A two component coupled wire model, viewed as 
a bilayer. (b) Schematic of tunneling processes in < |3.3|3.4|3.5[ ) 
that lead to the Moore Read state for bosons at filling v = 1. 
(b) After the transformations ( |3.6|3.7|3.12 1 the model de- 
scribes coupled strips of v = 1 Moore Read state coupled by 
tunneling electrons between the edge states, which are char- 
acterized by a c = 1 chiral charge mode and a c = 1/2 chiral 
Majorana fermion mode. 



The latter fields satisfy 



(3.8) 



for /i = p, a. 

For simplicity, we will first focus on the case in which 
tab = t, independent of a and b, and t, u and v are real. 
We will comment on the more general case later. In this 
case, in terms of the new variables, we have 



ah 



it a bOj t ab + h - c - = 

8tcos(4>f jP - #+i^)cos^cos$' +1)<7 .(3.9) 



and 



uO- 



h.c. = 2m cos 

IR 



+ 



3^ 



sin < 



2zw(sin 
h.c. = 2v cos( 
2u!(sin<^ cr sin^j^ 



— cos 



i^cos^J(3.10) 
^cos«^J(3.11) 



Note that in passing between the first and second lines of 
( 3.9 )-( 3.11) getting the factors of i right requires care in 



splitting the exponential. This is explained in Appendix 
A. 3, where the zero momentum components of 4>fi L are 
properly taken into account. 



Since the interaction term V is a sum of non commut- 
ing terms, analysis of this state is more complicated than 
it was for the Abelian quantum Hall states. However, a 
tremendous simplification occurs when the forward scat- 



tering interactions in are such that 

are decoupled, and the Hamiltonian for 



- h R/L 



and 



-R/L 
3," 



-R/L 
3," 

has the non 



interacting form, "W„ onint . in (2.7). In this case, the oper- 
ator exp has precisely the form of a bosonized Dirac 
fermion. This allows us to fermionize, by writing 



R/L R/L 



. R/L 



7H/t 
10 ■ 

e 



(3.12) 



where ip a is a Dirac fermion operator, and £j a , rjj a are 
Majorana fermion operators. For the charge sector, we 
define 



0j+l/2,p = {(t>f,p 

— (1R 



" + i,p)/2 

L 



<f>j+l/%a = (4>f,p + ^i+l,p)/2- 

They satisfy [d x 6^ p (x), 6 t > tP {x')\ 

[d x <fie,p{x),if£', p (x')] = and 

ld x 9e,p(x),ip e >,p(x')} = iir6i£i5 xx >. 
The Hamiltonian may now be written, 



(3.13) 



(3.14) 



n = n° Shh [d i!P ,ipi,p] + hUv*] + + v (3.15) 



where %g L L [@j,p 1 4>j,p\ has the form (2.8). For a special 
value of the interactions in the charge sector it is also 
possible to fermionize 9j p and <Pj p , though that is not 
necessary for our purposes. The free fermion Hamilto- 



R/L 



is 



nian for the Majorana fermion 77 ■ a 

HStfo*] =Y,j dxi (jfi&rfi* ~ Vjl&ri.) , (3.16) 

3 

with a similar expression for T-Lo\^j,a\- The interaction 
term is 



V = J2 J dx 



(u - v)i$At + (5 + v)iv**Vj,*\ ■ (3-17) 



We now assume 'Hg- L - L [9g p ,ip£ p \ 



is in a regime such 
that t is relevant, and 9e lP is pinned in a self con- 
sistent minimum of cos#£, p . % then describes inde- 
pendent free fermion problems for £ J:(T and T\^ a . The 
r\ sector has a gap with a k y independent dispersion 
E = ±^v 2 k 2 + (u + v) 2 . The £ sector has dispersion 

E = ±\/v 2 k 2 + \te ik y + u — v\ 2 with a gap that closes at 
a point when t = ±|u — v\ signaling a quantum phase 
transition. The phase diagram is shown in Fig. [5]i. 

We identify the t > \u — v\ phase with the Moore Read 
state^. Its physics is most transparent at the special 
point u = v where the chiral Majorana modes and 
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(a) MR : c = 3/2 (b) SP : c = 1 (c) B : c = 2 
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uV 



Tla Vi <t>p) uV 



4, 



> ) 



(d) t ab = t 




(e) t ab =tS ab 




FIG. 5. The upper panels are schematic diagrams depicting 
the coupling of the chiral edge modes. The solid lines rep- 
resent the charge modes <f>fp L - The dashed lines represent 

the chiral Majorana modes (f^ L and vfi L - Solid (dashed) 
arcs represent stronger (weaker) coupling, (a) describes the 
Moore Read (MR) state, with c—3/2 edge states, (b) de- 
scribes a strong pairing (SP) state with c = 1 edge states, 
and (c) describes a bilayer (B) state with c = 2 edge states. 
The lower panels are ternary phase diagrams as a function of 
t, u and v in the cases where (d) t a b = t, independent of a and 
b and (e) tab = t8 a b- The dashed line in (d) is the solvable 
line, where the decoupling of the chiral modes is perfect. 



ifi, L pair up with a pattern shown in Fig. In this 
case, it is clear that it has a left over gapless cmral charge 
mode cjJy p and a single chiral Majorana mode £,± a . This 
is precisely the structure of the edge of the v = 1 bosonic 
Moore Read stated Similar to the Abelian case, we may 
view the change of variables ( 3.6|3.7|3.12 ) as transform- 
ing a SLL of coupled bosonic wires to a SLL of strips of 
v = 1 Moore Read states coupled by their edge states. 
In this case, however, the coupling of the edge states 
goes in "both directions": the T^ a rfc a couples the edge 
states on a single strip, leaving behind the gapless edge 
states, while the cos 9j + i/2,p£f CT £j+i CT term couples the 
edge states on neighboring strips, leading to the bulk 
Moore Read state. 

As in the case of the Abelian quantum Hall states, 
bulk quasiparticle excitations are associated with kinks 
in 9t lP (x). The present case is slightly different, though 
because the transformation ipj i — > tpj.i + 27r (which con- 
nects equivalent states) translates to <l>jji a < ^p/cr + 7r - 
It then follows that the transformation 9j„ — > 9j tP + Tr/2, 
(£j,o-> Vj,cr) ~~ * ~(£j,<ri r lj,o-) connects equivalent states. 
The elementary quasiparticle is thus associated with a 
kink in which 9e tP (x) jumps by 7r/2, corresponding to a 
charge e/2. This introduces a domain wall where the 
mass term coupling and a changes sign. This 



binds a zero energy Majorana bound state, as is expected 
for the charge e/2 quasiparticles of the bosonic Moore 
Read state. As in Section II. C. 2, the quasiparticle tun- 
neling operators can be related to the backscattering of 
bare electrons. We defer the discussion of this to section 
IV.C. 

The t < |€t — v| phase corresponds to a strongly paired 
quantum Hall state of charge 2e bosons at filling v = 1/4. 
It is most easily understood in the limit t << \u — v\. In 
this case, the Majorana modes pair up with the pattern in 
Fig. [5]d, so that there are no gapless Majorana modes at 
the edge. In this limit, individual bosons can not tunnel 
between wires because it excites the gapped £j CT modes. 
However, a pair of bosons can tunnel without disturbing 
the £j ](7 sector. The charge modes thus pair up leaving a 
single gapless chiral charge mode at the edge. 

We now briefly discuss the case in which we relax our 
assumption about the equality of the different t a b- In this 
case, the t a b terms will have the general structure 



ab 



■h.c. = i ( if a rj 



L 



(3.18) 



where the 2x2 matrix T nm — t nm cos(20£. p +/3 a h) is char- 
acterized by a magnitude t nm and a phase (3 nm , which 
depend on i a6 . When 9 p is stiff, we may again analyze 
the problem by putting 9 p in a self consistent minimum 
and solving the non interacting fermion problem. 

It is clear that due to the existence of a bulk energy 
gap, the phases discussed above will persist for a finite 
range in the more general parameter space. However, 
another possible phase is possible where the neutral Ma- 
jorana modes pair up in the manner shown in FigJHJ;. In 
this case the edge has a gapless charge mode and two gap- 
less Majorana modes, or equivalently two gapless bosonic 
modes. 

To see this, consider another special limit where t a b — 
t5 a b, with t real. It then follows that t nm = t5 nm cos 26t lP , 
so that the V has a term fcos 2^ iP (^j ? (T ^ +lj(T + 
vfaVf+i a)- The phase diagram for this case is shown 
in Fig [5^. When u = v, the £ sector is gapped by t, while 
the 77 sector involves competition between t and u+v. For 
t < u + v we have the pairing in Fig. [H^i, giving the Moore 
Read state, while for t > u + v, we have the pairing in 
Fig. [5J:, which has two bosonic edge modes. This is most 
easily understood when t a f, — t8 a b and u = v = 0, which 
simply corresponds to two decoupled v — 1/2 bosonic 
quantum Hall states. 



B. Coset Construction 

The results of the previous section can be understood 
within the frame work o f the coset construction in con- 
formal field theory^^. This is useful because it helps 
us make contact with the work of Fradkin, Nayak and 
Shoutens^l as well as subsequent generalizations. It also 
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introduces a framework that will allow us to generalize 
our construction to the Read RezayP^ sequence of quan- 
tum Hall states. Here we give a brief introduction to 
this well developed, but somewhat abstract, mathemati- 
cal construction that emphasizes its physical meaning in 
the context of our coupled wire theory. 

Our construction began with the SLL fixed point, 
which has bosonic modes <pf^ L on each wire. The two 
right moving chiral modes on each wire define a confor- 
mal field theory with central charge c — 2. Through the 
sequence of transformations, the Hilbert space of these 
two modes was split into three pieces, described by 4>f pl 
and r\j, a - The bosonic mode corresponds to c = 1, 
while each of the Majorana modes has c = 1/2. At the 
point u — v > and t a b = t, the decoupling is perfect. 
The decomposition of the Hilbert space can be summa- 
rized by 



2 = 1/2 + 1/2 + 1. 



(3.19) 



The coupling terms we introduced allow the modes in 
the different sectors to pair up in different directions, as 
shown in Fig. [6] This leads to a bulk gap, but leaves 
behind edge states in some of the sectors. The Moore 
Read state thus has c = 3/2 at the edge. 

To understand this decomposition more generally, it is 
important to realize that at the special filling factor v — 
1/2 the original chiral modes (f>^ have an extra symmetry 
because expi<fi R has scaling dimension A = 1. It follows 
that the operators 



J* ± = e ± ^/(2nx c ), 



J 



(3.20) 
(3.21) 



generate an SU{2) symmetry. Each channel is thus de- 
scribed by a SU (2)i Wess Zumino Witten model. The 
two channels together have SU(2)i x SU(2)i symmetry. 
SU(2)i x SU(2) 1 has a diagonal subgroup, SU(2) 2 , gen- 
crated by J R = Jf + J.f. In terms of the boson operators 
we have 

jR = e ±< CQS ^R/^wXc), (3.22) 

J? = 9x</(2tt). (3.23) 

SU (2)2, in turn, has a subgroup [7(1) generated by Jf. 

The coset construction allows a Wess Zumino Witten 
(WZW) model described by a group G with subgroup 
H to be divided into two pieces described by G/H and 
H. This means that the Hamiltonian can be written as 
the sum of two commuting terms, Hg = Wg/h + Hh, 
so that the Hilbert space of eigenstates factorizes. In the 
language of conformal field theory the energy momen- 
tum tensor can be written Tq = T G / H + Th , and the 
components Tq/h and Th have no singularities in their 
operator product expansion. It follows that the central 
charge of a coset theory given simply by c G i H = cq — ch- 
Applied to the present problem, we have 

(3.24) 



SU(2) 1 xSU(2) 1 




SU(2) 1 xSU(2) 1 /SU(2) 2 
— gap 



SU(2) 2 /U(1) & U(1) 

< gap 



SU(2) 2 /U(1) & U(1) 
< — edge state 



FIG. 6. Coupling of edge states at the decoupled point 
described by the dashed line in Fig. [5ji. The right and 
left moving (7(1) charge modes and the 517(2)2/17(1) Ma- 
jorana fermion (Z2 parafermion) modes are coupled on neigh- 
boring wires, while the right and left moving SU(2)i x 
SU (2)i /SU (2)2 modes are coupled on the same wire. This 
pattern leaves U(l) and SU (2)2 /U(l) chiral modes at the 
edge. This provides a concrete interpretation for the coset 
construction for the Moore Read state. 



Using the fact that csu{2) k = 3fc/(fc + 2) it is simple to 
see that (3.24) is equivalent to ( |3.19 1. 

Consider the hopping term between wires, which for 
tab = t can be written 



ab 



0),ab 



— 8tJ + jJ_j +1 . 



(3.25) 



TsU(2) 1 xSU(2) 1 /SU(2) 2 + TsU(2) 2 /U(l) 



This term acts only in the SU{2)2 sector for the pair of 
chiral modes, and leads to an energy gap in that sector. 
Thus, if that were the only term (i.e. u = v = 0), then the 
SU(2) 2 /U(1) and [7(1) sectors would both be gapped, 
but each wire would retain c = 1/2 chiral modes from 
the S*C7(2)! x SU(2) 1 /SU(2) 2 sector. 

On the other hand, the term (u + v) sin^j^ sinc^^ 
acts only in the SU(2) 1 x S , C7(2)i/S , C7(2) 2 sector. Thus, 
if this is the only term (i.e. t — u — v = 0), then the 
SU(2) 1 x SU(2) 1 /SU(2) 2 is gapped, while there will be 
c = 3/2 chiral modes associated with the SU(2) 2 /U(1) 
and U(l) sectors. While this fact is clear from fermion 
representation ( |3.12[ ) , we will defer the proof that sin <p^ a 
acts only in the SU{2) X x SU{2) 1 / SU{2) 2 sector to sec- 
tion IV. B, where it will be demonstrated in a more gen- 
eral context. 

To summarize, the coset construction provides a way 
to fractionalizc a c = 2 Luttinger liquid into non trivial 
pieces. When those pieces pair up in "different direc- 
tions", as in Fig. [6j the resulting fully gapped phase is 
a quantum Hall state with edge states that reflect the 
non trivial coset conformal field theory. In the following 
section we will generalize this to develop a coupled wire 
construction for the Read Rezayi sequence of quantum 
Hall states, described by a SU{2)u theory. 

Before doing that, though, we will conclude this sec- 
tion by generalizing the coupled wire construction of the 
bosonic Moore Read state to describe fermions. 
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C. Generalization to fermions 

We now consider coupled wires of fermions. Unfortu- 
nately, for a uniform magnetic field there is no simple cou- 
pled wire construction for the v = 1/2 fermionic Moore 
Read state. The tunneling terms that are allowed by 
momentum conservation either lead to a strongly paired 
Abelian quantum Hall state of charge 2e bosons, or they 
involve pairs of non commuting terms that can not be 
easily analyzed using the present methods. Evidently, 
the Moore Read state is not sufficiently "close" to the 
SLL fixed point for uniform field. 

However, we found that if the magnetic field is stag- 
gered, so that the flux between neighboring wires alter- 
nates between two values, then a construction similar to 
the preceding section can be developed. One can view 
this as a generalization of the two channel construction 
in the preceding section, where instead of having the two 
layers directly on top of one another, one layer is slid over 
relative to the other. Equivalently, this can be viewed as 
a single layer system with a staggered field as in Fig. [7^,) . 
The system in Fig. [4^, has a magnetic flux per unit length 
in units of the flux quantum b = eaB/h that alternates 
between b and 0. Our more general construction then cor- 
responds to sliding one layer relative to the other, so that 
the flux per length in units of the flux quantum alternates 
between two values b% = b + 8b and b 2 = b — 5b. The aver- 
age flux b is related to the filling factor, v = 2kp/b. The 
two channel bosonic problem then corresponds to 8b = 6, 
while the uniform field corresponds to 5b — 0. We will 
show that when 8b = 2kp the allowed tunneling terms 
have a structure similar to that in ( 3.3|3.5 ). This con- 
struction gives the v — 1/2 fermionic Moore Read state 
for b = 4kF, as well as the more general "g-Pfaffian" 
state at v = 1/(1 + 17) for & = 2kp(l + q), where q is even 
(odd) for bosons (fermions) . The state in this series with 
q = — 1 is special, and corresponds to a p + ip supercon- 
ductor in zero net magnetic field. In our construction, 
the modification of the state by changing the uniform 
component of the field b is reminiscent of modifying the 
Moore Read wavefunction by including a Jastrow factor 
that compensates the change in magnetic fielcP. 



1. q-pfaffian state 

Consider an array of wires with alternating magnetic 
flux, shown in Fig [7] We parameterize the two fluxes 
as 2kp(2 + q) and 2k F q. We group the wires into pairs, 
indexed by j and a = 1,2. The interaction terms then 
have a form similar to ( 3.3|3.5 ): 



j J \ab=l ) 

(3.26) 

There are four terms coupling pair j to j + 1, 

Qt ^ _ e i(^j,a-Vj + l,b) + (q+2)(,e j]a +6 j + 1]b ))Q^ ^ (3.27) 




(b) 



j.1 



0/ 



1.2 



j+1,1 j+1,2 



2q 



c:> 



cn> 



2q 



c:> 



cb 



C^> 



c:> c:> 



c:> 



cp q 



C^> c~> 



2q 



2q 



2 +q yOj +1/2 ,ab 



CD 



FIG. 7. (a) Coupled wire construction with a staggered mag- 
netic field. For Sb = b, the theory is equivalent to the two 
component boson system shown in Fig [4^. (b) Schematic di- 
agram showing the coupling of the edge states, similar to Fig 



with 



1 e i2 <3 e j + i,i 



(3.28) 



Two terms operate within a single pair. The first involv- 
ing tunneling an electron between the two wires 



o i(<Pj,i-<Pf,z+q(6j,i+0j,a)) 



(3.29) 



The second giving an interaction between the 2kp densi- 
ties. 



qv _ e i(%,i-%,2) 



(3.30) 



From (2.16) and (2.261, it is clear that these interac- 



tions are appropriate for bosons (fermions) when q is even 
(odd). 

We can write the first term as 



0* , = e %,a-*i+l,t 



(3.31) 



with 



ft 
ft 
ft 
ft 



(p jtl + (q + 2)9 hl + 2q9 jl2 

~ il + 2 )0j,i 
<p ji2 + (q + 2)^-2 
<Pi,2 - (l + 2)0j,2 - 2gfl 3 -,i 
We next define the sum and difference variables 



(3.32) 



~ h RIL 

l R ' L 
3,a 



Iff 



ft L )/2, 
lfi L )/2. 



(3.33) 
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These variables obey the commutation relations, 



d x 4>l p (x),4> p f Jx') = 2m(l + q)p8p P '5 3J ,8 xx {Z.M) 



d x $ (x), $j,Jx') = 2mp6pp'5 jj <6 x 



Note that the commutation relation for the a sector is 
identical to (3.8 1. This allows us to proceed in the same 



manner as section III. A. For the a sector to be unaltered, 
it was essential that the staggered field satisfy Sb = 2k p. 
The charge sector, on the other hand is modified, and 
resembles that of the Laughlin state with m = 1 + q. 

As in section III. A we focus on the case where t a b = t, 
independent of a and b. Then 



E 



^W C OS^COS^ + l,a, 



along with 



(3.36) 



(3.37) 
(3.38) 



From this point the analysis proceeds exactly as in sec- 
tion III. A and will not be repeated. The only difference is 
that the charge sector has a modified structure constant 
in (3.34 1. This changes the exponents for tunneling elec- 



trons into the edge states, as well as the charge of the 
bulk quasiparticles. As expected, for the v — 1/2 Moore 
Read state (q = 1), the quasiparticles have charge e/4. 



2. p+ip superconductor 

The special case q = —1 in the previous section cor- 
responds to a fermion system with 6 = 0. This is not a 
quantum Hall state, but rather a. p + ip superconductor^, 
and the analysis is slightly different. For q = — 1, the 



transformation in (|3.32|) breaks down because <f>^i% 



come linearly dependent. This can be seen from the fact 



that <fij p and 4>j p , defined in (3.331 are in fact the same 
operator. 



We proceed by defining cj)^ and (j>f a as in (3.321, but 



replacing tpf and <j& by 



<Pj,P 



0j,i -+ 



- <Pj,2 



+ 9(^,2 



))/2 



(3.39) 



These satisfy [d x 9j. p (x), ipj> p (x')] = in5jj>5 xx i. It follows 



that 4>"[p L = <fij, P ± (1 + <l)6j,p- For q = — 1 the analysis 
is then the same, except that ( |3.36 1 becomes 



E °S,<* = cos^ cos 

ab 



bj + l,a- 



(3.40) 



This term has precisely the form of the tunneling of 
electrons between the edge states of strips of p + ip su- 
perconductor. When O* ab flows to strong coupling the a 



Josephson coupling cos2(0j- jCT — <f>j + x,a) wm be generated, 
and the phases ipj iP on neighboring will lock together. 
However, unlike the quantum Hall case, there will be a 
gapless bulk collective mode associated with slowly vary- 



(3.35) ing fluctuations in tpj tP - The neutral sector is identical 
to that of the Moore Read state, however, for the charge 
mode, the gapless edge mode in the quantum Hall case 
is replaced by a gapless bulk collective mode. 



IV. READ REZAYI SEQUENCE 

In this section, we will generalize the coupled wire 
construction to describe the Read Rezayi sequence of 
statesP. This sequence includes the Moore Read state 
for k = 2, as well as other states, which are described 
in terms of the parafermion conformal field theory^. 
As in the previous section, the analysis is simplest for 
bosons. Following the analysis of Fradkin, Nayak and 
ShoutemP, we thus consider k channels of bosons, which 
are each at filling v — 1/2, so that the total filling factor 
is k/2. At the end of this section we will briefly describe 
the generalization, similar to section III.C, which gives 
the known Read Rezayi states at v — k/(2 + kq), where 
q is even (odd) for bosons (fermions). 



A. Bosons at v = k/2 

Consider coupled wires of k channel bosons. The 
analysis is similar to Section III. A, except now each 



wire is characterized by 9j, a and (pj <a satisfying (3.1) for 
a = l,...,fc. The Hamiltonian is again 'Hg LL (0, <j>) + V 
with 

(k 
E + 
ab=l 

k \ 

]T u ab Ol ab + v ab Ol ab ) +h.c. 



a<b=l / 

The interaction coupling neighboring wires 

H* — p i (.'Pi,a-Vi+l,b+2(9j,a+8j + l,b)) 



(4.1) 



(4.2) 



is the same as before, while the interactions operating 
within a single wire come in more varieties, 



U j,ab 



oK<Pj,a-<flj,b) 



(4.3) 



as well as an interaction that locks the "2fcp>" densities 
of the two channels, 



U j,ab 



e i(20j, o -20 3 ',i,)_ (44) 

As in section III. A we first define chiral boson modes 

4>? <a = (fj,a + 20j,a)/V2 

4>la - (<Pj,a - 20 J -, a )/V2. (4.5) 
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3.6 



For later convenience, this definit ion d iffers by a factor of 
\J 2 from the modes defined in Eq. 
a charge 
writing 



We then introduce 

a charge mode 4>fp L and k — 1 neutral modes $f^ L by 



l R ' L = 

3-M 



~R/L 

V 3,P 
7R/L 



with 



-R/L 
V 3,V 



R/L 



Ofj, a is an orthogonal matrix that has the form 



(4.6) 



(4.7) 



(4.8) 



where d a are a set of k vectors with k — 1 components 
that satisfy 



a 

£« = <w 



1/fc. 



(4.9) 
(4.10) 
(4.11) 



d a may be viewed as the unit vector in the a direction 
projected into the plane perpendicular to (1, 1, 1). For 
example, for k = 3 they form 3 planar vectors oriented 
at 120°, 








</, = ( -| ) .(4.12) 



The transformation then has the explicit form, 

k 



~ b R,L 
3,P 



if- 



R/L 



a=l 



along with 



3,a 



1 ~ 



R/L 



(4.13) 



(4.14) 



For k = 2 the charge and spin modes <pf^ L are identical 
to the corresponding modes defined in Eq. |3.7|in Section 



III. A (though 4>^ a L differ by y/2. The charge and neutral 
modes satisfy 

d x4>^ fl (x),^ /tfl ,(x') =2ivipSp P >8^ li 'Sjji6 xx > (4.15) 

Expressed in these variables, the interaction terms 
have the form, 



n v _ j(i a -dt b y($f :tr -^)/V2 

^j,ab ~ e 



(4.17) 
(4.18) 



We will now focus of the special case in which t a b = t 
are independent of a and b. Then, 



3 

£K& + v ab )iTf ab T^ ab + (u ab - v ab )iZf tab Zfy.l9) 



ab 

where 



and 



iV2d a ■ 



r^R/L 

T 3,ab = Sm 



~.R/L 



~ d b ) ■ $f/ a L 



"j.ab 



- (d a - d b ) ■ (fi 



R/L 



(4.20) 

(4.21) 
(4.22) 



In the following we will show that at the special point 
u = v a decoupling similar to what occurred in Section 
III.B occurs. To establish this, we will first use the coset 
construction to show how the k chiral modes on each wire 
decouple into separate sectors. We will then show that 
^f^ L acts only in one sector, while acts only in the 



other. The coupling terms in (4.1 1 then lead to gaps in 



which the different sectors are paired in different direc- 

Rl L 

tions, leaving behind non trivial edge states. \P ■ will 
be identified as a Z^ parafermion operator. The coupling 
term Y^ ah T^ ob , on the other hand, leads to a theory on 
an individual wire which can be identified with the crit- 
ical point of a Z-^ model, which is a particular k state 
generalization of the Ising and 3 state Potts model. 



B. Coset Construction and Primary Fields 

Each wire is characterized by k right and left mov- 
ing chiral modes, which individually are described by 
a SU(2)i WZW model. As in the previous section, 
[SU (2)i} k can be decoupled by considering the diagonal 
sub algebra SU(2)k- This leads to the following decom- 
position of the energy momentum tensor 



T, 



[SU{2)i] k 



(4.23) 



T [SU(2) 1 ]"/SU(2) k + T S U(2) k /U(l) + T U{1) . 

In terms of the central charge, this is equivalent to 



k = 



k{k-l) 2(k-l) 
k+2 + k+2 



1 



(4.24) 



Clearly k = 2 reduces to (3.191. For k = 3, we have 
3 = 6/5+4/5+1. The SU{2y k ~JU(l) sector is precisely the 
Zk parafermion theory introduced by Zamolodchikov and 
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[SU(2) 1 ] k /SU(2) k 
< gap 




SU(2) k /U(1) 

< gap 



& u(i; 



SU(2) k /U(1) & U(1) 
I — edge state 



FIG. 8. Along the solvable line u a t = v a b the right and 
left moving (7(1) charge modes and the SU(2) k /U(l) Z k 
parafermion modes are coupled only on neighboring wires, 
while the right and left moving [SU(2) 1 } k /SU(2) k modes are 
coupled only on the same wire. This pattern leaves U (1) and 
SU(2)k/U(l) chiral modes at the edge. This provides a con- 
crete interpretation for the coset construction for the Read 
Rezayi state. 



Fateev^. The c = k theory thus decomposes into a [7(1) 
charge sector, a SU(2)k/U(l) parafermion sector and a 
third sector described by the coset [SU(2)] k / SU{2) k . 

To show this decoupling explicitly, and to demonstrate 
that ^ and T are primary fields that act only in the de- 
coupled sectors, we explicitly construct the energy mo- 
mentum tensors. Here we consider the right moving chi- 
ral sector on a single wire j and omit the superscript R 
and subscript j. The chiral channel is described by the 
Hamiltonian 



H= dx 



'47T 



((d z p ) 2 + (d z $ a f) . (4.25) 



In order make contact with the conformal field theory lit- 
erature, we focus on the energy momentum tensor, which 
is closely related to the Hamiltonian density, 



1 



(4.26) 



The translation between the Hamiltonian formulation 
and the conformal field theory is briefly reviewed in Ap- 
pendix B, where we also show t hat T \ SU ( 2 )^k decomposes 
into three terms, according to (4.23), with 



r. 



SU(2) k /U(l) 



k + 2 



(4.27) 



I (4-28) 



and 



'[SU(2) 1 ]"/SU{2) k 



k + 2 \ 2 



2 e iV2(d a -d b y 



.(4.29) 



The non trivial content of this decoupling is that the 
three energy momentum tensors have no singular terms 
in their operator product expansion. In a Hamiltonian 



formalism, this means that the Hamiltonian (4.25), splits 
into three commuting pieces, H = Hu(i) +T~Lsu{2) k /u{i) + 
T~L[su(2)!\ k /su{2) k ■ This decoupling has also appeared in 
a somewhat different context in Ref. [55] 

We now show that if? and T act only in a single sec- 
tor. This is done by computing the operator product 
expansion with the T's. Details of the calculation are in 
Appendix C. We find that for z — > w the singular terms 
in the operator product expansion are 



T S u(2) k /u(i)(z)^(w) 
1 - l/k 



(z — 11))'* 



1 



(z — w) 



d,*(z), (4.30) 



T u(1) {z)*(w) = T [su(2h]k/su{2h ^(w) = 0. (4.31) 



Eqs. |4.30| and |4.31| shows that <f is the primary field of 
the SU(2) k /U(l) theory with scaling dimension 1 — l/k, 
known as a Z k parafermion operator, f is a gener- 
alization of the Majorana fermion, which can be re- 
garded as a Z 2 parafermion. The fact that there no 
singular terms in the OPE for the other two sectors 
means that if? acts only in the SU(2)k/U(l) parafermion 
sector. In a Hamiltonian formulation, we would have 
= [*,^[si/(2)i]*/sri7(2) fc ] = 0. A mass term 
vJr fl tij/£ h as scaling dimension 2 — 2/fc, and is relevant. It 
leads to an energy gap in the parafermion sector. 
For T a b we find 



T [SU{2) 1 ] k /SU(2) k (^)Tab(w) 

1/2 



(4.32) 



(z — w)' 



rT afc (z) + 



(z — w) 



d z T ab (z), 



T, 



u(i)(z)Y ab (w) = T su{2)k/u{l) {z)T ab (w) = 0. (4.33) 

This shows that T ao are primary fields of the 
[SU{2)} k /SU{2) k sector with scaling dimension 1/2, and 
do not act in the SU(2) k /U(l) or the 17(1) sectors. A 
mass term T R T L has dimension 1 and leads to an energy 
gap in the [SU(2)] k / SU(2) k sector. 

We have also computed the OPE's for S &, defined in 
(4.22). Unlike ^ and T, though, 3 is not primary and 
acts non trivially in both the [SU(2)] k / SU(2) k and the 
SU(2) k /U(l) sectors. Thus, unlike the k — 2 case, it 
is not clear how S fl S L competes with the other terms. 
Nonetheless, on the special line u ab = v ab , t ab = t the S a f, 
term is absent, and we have the decoupling shown in Fig. 
|9j in which the U{1) and SU(2) k /U(l) sectors are gapped 
across wires, while the [SU(2)] k / SU(2) k sector is gapped 
within a wire. This gives the Read Rezayi state, which 
has a left over gapless edge state with a U(l) charge mode 
and a SU(2) k /U(l) Z k parafermion mode. 
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C. Quasiparticle Operators 

To construct quasiparticle operators we follow the logic 
of Section II. C. 2 and consider the 2fcp backscattering of 
bare particles on channel a of wire j, 



Xj,a 



„2i6j 



(4.34) 



= eVW^LH^-^-^) (4.35) 
We thus define 



yR/L = R/L 



where 



S ,,a = ex P 



iV2d a ■ 



TR/L 



(4.36) 



(4.37) 



R/L 



QPj+i/2a crea tes a quasiparticle with charge 1/2 
(check!) with a non trivial action in the neutral sector. 

Like S a f,, the operators S a are not primary, and acts 
in both the [SU(2)} k /SU(2) k and SU(2) k /U(l) sectors. 
Nonetheless, in the next section we will argue that when 
the [SU{2)] k / SU(2) k sector is gapped, E a acts as a pri- 
mary field projected into the parafermion sector, which 
corresponds to the spin operator a. 



D. Relation to Z k Statistical Mechanics Model 

On a single wire, the mass term (u ab + v ab )T^ b T^ b 
opens a gap and leaves behind a SU(2) k gapless edge 
state with a charge mode and a Z k parafermion. The 
Z k parafermion conformal field theory is known to de- 
scribe the critical point of a Z k generalization of the Ising 
modeP^l. For k — 3 it is the 3 state Potts model. Simi- 
larly, for A; > 3 it is a particular version of a Z k symmetric 
k state generalized Potts model. In this section we show 
that our bosonized representation provides a simple and 
intuitively appealing way to understand this connection. 
This allows us to identify the projected operators S a with 
the primary fields a of the Z k parafermion model. 

We begin by rewriting the mass term by introducing 
new variables for the single wire, 



These variables satisfy 



#)/2. 



[d x eS(x),<4tf)]=in5 aP 6 x 



(4.38) 



(4.39) 



The Hamiltonian for the neutral sector of a single wire 
then has the form 



H 



Z7T V 

^ u ab cos V2d ab ■ ^ + v ab cos V2d ab ■ $4.40) 



O 
A 



O 
B 



O 
A 




O 
A 



O 

c 



o 

A 



FIG. 9. Lattice of minima of the periodic potential 
— ~}2 ab cos \ /r 2d a t ■ 6 a as a function of 8 a . The shaded region 
shows the set of distinct values of 6* CT , which is compactified to 
a torus due to the vortices created by expi\> r 2d a t > ■ <p a . There 
are thus three distinct, but equivalent minima, labeled A, B 
and C. 



where we have abbreviated d ab = d a — d b . 

First consider the simplest case k = 2, where 9 a and 
ip a have a single component, and d\ 2 — d 2 i — v2- Then 
we have 



H 



2tt 



+ {d x ip a ) ) — u cos 28 a — v cos 2tp a 

(4-41) 

Viewed as a transfer matrix for the partition function of 
an anisotropic statistical mechanics problem, this Hamil- 
tonian gives a well known representation of the 2D Ising 
modePSl. This can be understood by first considering 
u = 0. This describes the 2D XY model with order pa- 
rameter (cos 9 a , sin^ CT ). For v = 0, 9 a is a non compact 
variable, so there are no vortices. From (4.39) it can 



a b 



be seen that exp ±2i(p a (x, r) creates a vortex where Re- 
winds by ±27r around (x,r). v is thus the fugacity for 
vortices, and its presence makes 9 a an angular variable 
defined modulo 2ir. Integrating out 9 a gives the sine 
gordon representation of the XY model. For nonzero 
u, — cos 29 a introduces an Ising anisotropy into the XY 
model. For large u 9 a is pinned in the minima of this 
potential at 9 = nn. Due to the presence of v, only two 
of those minima are distinct. For u ^ v, since both u 
and v are relevant, the system flows at low energy to a 
strong coupling phase in which either 9 a or <p a is pinned. 
The symmetry under u 4-> v and 9 a <-> <p a is precisely the 
Kramers Wannicr duality of the Ising model. At the self 
dual point u — v, the system at low energy flows to the 
fixed point of the Ising critical point. 

For k > 2, a similar interpretation is possible. Now, 
however, 9 a lives in k — 1 dimensions. expi\/2d ab ■ (p a 
creates vortices around which 9 — >• 6 a + ^/2Ttd ab . This 
compactifies 9 a , so that it is defined on a k — 1 dimen- 
sional torus, cos \f2d ab -9 a introduces a periodic potential 
for 9 a, and pins 9 in its minima. 

For k = 3 the minima of the periodic potential are 
shown in Fig. [9] They form a triangular lattice with 
lattice constant 27r/v3. The compactification of 6 a is 
associated with a larger triangular lattice with lattice 
constant 27r. This identifies points on the original lattice 
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that are on the same y/3 sublattice. There are 3 distinct 
minima. We thus have a two dimensional generalization 
of the XY model (where the order parameter is defined 
on a torus T 2 ), with a 3 state anisotropy. Since both 
the vortices and the anisotropy are relevant (the peri- 
odic potential has scaling dimension 1), this leads to a 
3 state generalization of the Ising model with Z 3 sym- 
metry, which is uniquely specified by the 3 state Potts 
model. Again, the critical point appears at the self dual 
point u ab = v ab . 

For k = 4, the minima of the periodic potential form 
a three dimensional fee lattice. The fee lattice can be 
viewed as a larger bee lattice with a 4 site basis. The 
compactification is associated with the larger bec lattice. 
There are thus 4 states, so we have a 4 state generaliza- 
tion of the Ising model with Z4 symmetry, which is known 
as the Ashkin Teller model. This model has a parameter, 
which for different values gives, for example, the 4 state 
Potts model and the 4 state clock model. It is not imme- 
diately obvious what the value of that parameter should 
be based on the form of (4.41 1. However, from the analy- 



sis of the previous section, we know that the critical point 
is described by the Z4 parafermion conformal field the- 
ory. We thus expect that this model describes the Fateev 
Zamolodchikov point of the Ashkin Teller model. 

For general k the minima occur on a k — 1 dimensional 
lattice formed by combinations of \/2-nd a . There are k 
distinct but equivalent minima to this potential, which 
can be located at 9 a = 9 n = nv^irdx, for n = 0, k— 1. 
The minima at n — k is equivalent to the one at n = 
because from (4.9) kd\ = Y] ^ d a j . All other minima 

of cos \[2d a b ■ Oct can also be reduced to these k minima 
with a suitable combination of \f2i:d a ) ] . This model thus 
describes a k state system with Z k symmetry. Z^ models 
have extra parameters for k > 4, but as discussed above, 
since the critical point is described by the Z^ parafermion 
theory we conclude that our model describes the Fateev 
Zamolodchikov point of the Z^ model. 

Now we consider the operators E a discussed above, 
which has a simple interpretation. When (u a f, + 
Vab^ab^ab opens a large gap, then we can restrict 9 a 
to the minima 9 n — nyhKd\. It is then straightforward 
to see that 
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FIG. 10. Schematic diagram for the generalized Read Rezayi 
state at filling v = k/(2 + kq). groups of k wires are coupled to 
one another by Oj +1 / 2 aby which require a specific staggered 
magnetic field for momentum conservation. Wires within a 
group are coupled by OJ'J 6 , which are independent of the field. 
Representative examples of 0* +1/ / 2 a b an d Ojab are shown. 



E. Generalization 

As in Section III.C our construction for the Read 
Rezayi sequence can be generalized by introducing a stag- 
gered component to the field. Again, the way to think 
about it is to start with the bosonic state at v = k/2, 
which can be viewed as a staggered field, in which the 
field between neighboring wires is b for one out of every 
k neighbors and for the other k — 1 neighbors. Keep- 
ing this staggered field fixed, we now add a uniform field 
6, and find that for certain values, which correspond to 
filling factor 



k 



kq 



(4.43) 



there are allowed tunneling processes, which have a 
structure similar to ( 4.1||4.4 ). Expressed in terms of 
charge/neutral variables, as in (4.13), we find that the 
neutral sector is independent of q, while the charge sec- 
tor is modified, as in (3.34). 



Rather than repeating the algebra in Section III.C, we 
simply display the diagram, analogous to Fig [7j in Fig 

M 



CONCLUSION 



iR yL _ i\p2d a -e n _ -2mn/k 
'j,a j,a ' 



independent of a. This is precisely the spin order pa- 
rameter of the Zfe model, which gives different values 
e -27rm/fc £ or different ordered states specified by n. 
We thus conclude that when the [SU{2){\ k / 'SU(2) k sec- 
tor is gapped, the operator Sj^ a , when projected into 
the SU(2)k/U(l) sector corresponds precisely to the spin 
field a of the Zp. model. 



In this paper we have introduced a new formulation 
of non-Abelian fractional quantum Hall states, in which 
electronic models built from coupled interacting one di- 
mensional wires can be analyzed using Abelian bosoniza- 
tion. The picture that emerges from this analysis is sum- 
marized in Figs. [6] and [8] Non-Abelian states can be 
viewed as systems in which the original one dimensional 
chiral fcrmion modes are split into fractionalized sectors, 
in accordance with the coset construction of conformal 
field theory. The different coset sectors are then coupled 
to one another in "opposite directions" . This leads to 
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a simple understanding of the edge states and quasipar- 
ticles that is similar in some respects to the one dimen- 
sional AKLT model. In the case of the Moore Read state, 
the c — 1/2 coset theory has a simple fermion representa- 
tion. For the more general Read Rezayi states, the coset 
theory can be identified by a mapping to the critical point 
of a Zk statistical mechanics problem. We now conclude 
with some future directions and open problems. 

It is natural to speculate that a coupled wire construc- 
tion is possible for any quantum Hall state. For example, 
we expect that it should be possible to construct the hi- 
erarchical generalizations of the Moore Read and Read 
Rezayi states^. In addition, it would be interesting to 
develop the construction for the non-Abelian spin singlet 
state introduced by Ardonne and Schoutend^, which is 
based on a 5/7(3)2 coset theory. More generally, a much 
wider variety of quantum Hall states can be formulated 
using the parton construction. It would be interesting 
to understand the connection between our more concrete 
approach - which can be formulated in terms of electrons 
and solved - and the more abstract constructions. Our 
approach should also be applicable to spin models, and 
it is likely an anisotropic version of Kitaev's honeycomb 
lattice modef^l could be analyzed. 

It would be interesting to further explore the way in 
which the coupled wire construction accounts for the non- 
Abelian statistics of the quasiparticles. In the Abelian 
case, as explained in Section II. C. 2, the Abelian statistics 
of the Laughlin quasiparticles can be simply understood. 
It should be possible to understand the degeneracy and 
braiding properties associated with quasiparticles in the 
non-Abelian case in terms of the coupled wire model. 

Finally, we ask whether the coupled wire construction 
could be useful for describing any phases which are not 
already well understood. It seems likely that the cou- 
pled wire model could give a concrete representation of 
fractional Chern insulators, non-Abclian Chern insula- 
tors and fractional topological insulators in two, and pos- 
sibly three dimensions. 
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Appendix A: Klein Factors and Zero Modes 

In this appendix we carefully treat the Klein factors 
for fermions, along with the zero momentum modes of (p 
and 8. This shows that in expressions such as (2.331 and 



be treated separately. Define 



(2.55), the Klein factors may be ignored. 



An explicit formulation of the Klein factors can be for- 
mulated by considering a system with periodic boundary 
conditions in the x direction. The chiral modes may then 



= 7= — e 
y/2nx. 



(Al) 



Here, p = R/L = +1/ — 1 describes the right and left 
moving chiral modes, and 



[4%(x),</%,(x )] = ipirSpp'Sjj, 



(A2) 



where s xx i — sgn(a; — x'). The Klein factors are neces- 
sary to assure that fermion fields associated with different 
channels anticommute. They may be represented by the 
factors 

^ = (-1)So„x«,») jv j 1 (A3) 
where the number operator for each chiral channel 

Nf =p [ dxd x tf/2n, (A4) 



satisfies [Nj,^,] = i5jj'5p P ' and has integer eigenvalues. 
We have chosen an ordering of the chiral modes with 
direction p = R/L on wire i, such that (i,L) < (i,R) < 
(i + 1,L) < (i + 1,R). Defined in this way, the Klein 
factors mutually commute [kP-, kP-,] — 0, but the fermion 
operators anticommute, {tpji^Pf} = 0- Other choices for 



the phase factors in ( A3 1 are also possible 



The density and phase fields defined on each wire may 
be introduced as 



^• = (^ + #+7^/0/2 
9j = (<f>f-tf+TrNh/2 



(A5) 
(A6) 

These satisfy [0j(x),9j>(x')] — \tpj(x),(pj>(x')] = 0, along 
with 



[9j(x), <Pj'{x )] = inSjjiQ(x — x) 
The electron operators are then 



\/2ttx c 



(A7) 
(A8) 



where the Klein factor (with no superscript ) used in ( 2 . 7 1 , 
is 



Ki = (-1) S *< 



(A9) 



is now independent of p. 

Consider the backscattering operator on an individual 



wire, 



2-KlXr 



(A10) 



The Klein factors are absent, and can safely be ignored. 

In the subsections 1 and 2 we will apply this anal- 
ysis to the Laughlin states and hierarchy states. The 
bosonic Moore Read state does not require Klein factors 
in the original model, however, they require care when 
refermionizing. This is discussed in subsection 3. 
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1. Laughlin States 



2. Hierarchy States 



The electron operator responsible for the fermionic 
Laughlin state v = 1/m, for m odd, is 



O 



X/r. 



We write this as 



1/m - ib L ^ \b 



where 



= W-f) 



■(#) 



(All) 
(A12) 

(A13) 



We now keep the Klein factors and define Oj+1/2 and 
<Pj+x/2 such that 



Then, 



where 
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2^5 



3+1/2 
'3 + 1/2 



^+1 

%+i 



(A14) 



(A15) 



1 + m 
2 

1 — TO 



1 — TO 
1 + TO 



(A16) 



are sums of N^ L determined by the Klein fac- 



tors, using ( Al|A3|Al3 ). Defined in this way, the 
commutation relations obeyed by ipf^ L guarantee that 
j e 2i^(x) )6 2is,,(x')] = for A,B = §,<p and x ^ x'. This 
means that [A e (x), B e > (x')] = iP^Pn/2, where P^P is 
an integer. However, there is freedom in how is dc- 



fined because (|A3j) is unchanged when N- is increased 
by an even integer (which could depend on N^ L ). This 
freedom can be exploited to define 9 and <j> so that they 
obey a standard commutation relation. While a general 
method for determining remains to be developed, 

we have by trial and error found (non unique) solutions. 
For 



N = 
3 



1 



~-N 



N 



1 



TO — 1 



the fields 6*,- 



2 

and 



2 

TO 



-N 



j+x 



-N 



j+X) 



(A17) 
(A18) 



defined in (A15) satisfy 
jAl|),([A3l), JAI3| ) and ( |A14| ), as well as [6i(x), 9j(x')] = 



wKx),(f)j{x)\ = and 

= wr£y8(a: - x'). 



(A19) 



The procedure for defining the Klein factors for the 
hierarchy states is similar to that in the preceding section. 
Here we just sketch the process. We may again write the 
tunneling operators, defined in (2.45) as 



02k = ^U,X^k 



with 



fa 

4*2 



(^-1)^(^-1)^(^2^21) 



W4-1 

(*4)"* 



i^lx 



(A20) 



(A21) 



i>k,2 = (^ t )^(^ T )^(«-i^ fe -i) 



We now define 



^ = rttxJL 



Then, 



with 

4>lx 



4>k,2 

and 



26 

2y5 



fe+1/2,0 
fe+1/2,0 



t>a,k+X 
^a,k+l 



n + mi R 

2 ™2fc-l 

n — mi R 



n-mi L 

2 ^2/c-l 

n + toi £ 



2 


V2fc- 


2 


Y2k-l 




n + mi 
2 




n — mi 
+ 2 < 


">2k 




n — toi 


& 


n + TOi 


*>2k ~ mo(<p2k 




2 




+ 2 < 


-1 


JV fe+l/2,l - 


n - 


" 2 W1 (^- 


-1 + ^2fc +1 ) " 






(n 


+ 


- + 




JV fe+l/2,2 ~ 


n - 


" 2 W1 (^ 








(n 




H - <^2k + 




iV fc+i/2a 


n - 




-1 - N£ k+1 ) - 






(n 




- " 




Jv fc+l/2,2 


n - 




- fl» +2 ) + 






(n 




- n(iV« - N? k+1 



(A22) 



(A23) 



4) 



(A24) 



02fc-l)- 



(A25) 



)• 

These fields then satisfy 

[0k.a{ x )^l,b( x ')} = [ ( ^fc,a( a; ) I W,b(^')] = °- 

In addition 

[e k , a x, <fi,b(x')} = iir5 kl (K ab e(x - x') + W ab ), (A27) 



(A26) 
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with K given in (2.491 and W ab are integers. It may be 



possible to find another choice for N k ' in which W ab = 
0. However, for the purpose of this paper, this choice 
suffices. 



3. Moore-Read state 

Here we provide the details of the zero modes and Klein 
factors for the v — 1/2 bosonic Moore Read state dis- 
cussed in Section III. A. Since it is a model of bosons, 
there are no Klein factors in the original model. How- 
ever, it is necessary to keep track of the zero modes in 
order to correctly fermionize the theory. The analysis 
in this appendix leads to the appro priate (an d slightly 
counterintuitive) factors of i in Eqs. ( 3.9|3.11 1. 

Recall $J a (x) ~ e lVi < a ^ is the boson creation opera- 
tor and p n ~ e m{2kfx+2e(x)) re p resen t s the density wave 
at q ~ 2'Kpn. In order for these fields to commute for 
i/i', we require 



[0j t a(x),(pf tb (x')] = iTrSjjiS ab Q(x - x') 



(A28) 



where Q(x — x') is the step function. To define the chiral 
modes, with the appropriate commutation relations, we 
must augment (3.6) with the appropriate factors as in 



( A6 ) . This may be written 



u J,a 
<Pj,a 



Here, N? a 



P Jdxd x 4>lJ(2n) 



- nN£ a )/2. (A29) 
These fields obey com- 



mutation relations identical to ( A2 1. The extra N^ a term 
was necessary to make 4>^ a and (j>^ a com mute . Charge 
and spin modes can then be defined, as in (3.7 1, and the 



correction terms involve Nf ,, 
satisfy 



p J dxd x (p p j fl /(2Tr), which 



(A30) 



for [i — p,a. 

In terms of these variables we then write 



= exp[i( 
O" = exp[i( 



' 3 ~~*-l-\-r],<T tj,<t ■ ■•'• 3 ,<r/l- ( A31 ) 

We now separate the chiral modes in the exponential, 
keeping track of the commutator, and obtain 



(A32) 



It then follows that 
0V + h.c. 



i[cos(4>f t0 . + ttN^) cos 4> 



3,° 

sin 

R 



3,° 



f„ + nN^) sin & 



O) + h.c. = i[cosO£ CT + ttN^) cos </> 



3, a 



(A33) 



• ttJV^J sin 



Thus, the Josephson and inter wire coupling terms are 
given in (3.1013.11 ). A similar analysis gives (3.9). The 



factor of in (|A31|) provides precisely the necessary 



factor to define the Klein factor for the fcrmions in (3.12 ), 



as in Eqs. AT and A"3 



A similar analysis can be applied to the Read Rezayi 
states for general k. For example, (4.17) is modified to 
be 



(A34) 



where N p a = pj dxd x (/> P a /2ir and d ab = d a - d b . 



Appendix B: Decoupling of Energy Momentum 
Tensor 



In this appendix we review how the energy mo- 
mentum tensor for [SU(2)i] k is decomposed into 
[SU{2) 1 } k /SU(2) k , SU{2) k /U{\) and U{1). 

Since the notations in the conformal field theory liter- 
ature and the condensed matter literature are somewhat 
different, we first review the translation between the two 
for a single non interacting (right moving) fermion mode. 
In the CFT notation, this is written as 



(Bl) 



where z is a complex coordinate and the colons denote 
normal ordering. The chiral boson field <p satisfies 



(cj){z)(j)(w)} = - log(> - w) 



(B2) 



so that if> satisfies (ijr (z)ifj(w)) — l/(z — w). The energy 
momentum tensor is 



T(z) 



1 



(B3) 



Using (B2) and Wick's theorem, it can then be shown 



that the singular terms in the operator product expansion 
(OPE) of T(z) with ip(w) are 



T(z)ip(w) 



(z — w) [ 



1 



-d w iP(w). (B4) 



with A = 1/2. This shows that if) is a primary field with 
dimension A = 1/2. 

In the condensed matter literature the chiral fermion 
field is often written as 



ip(x,T) 



1 

\J2itx c 



J<j>(x,T) 



(B5) 



where the operator is not normal ordered and x c is 
a convergence factor in divergent momentum integrals, 
which plays the role of a short distance cutoff. Since 
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expi0 = (xcjL) 1 ! 2 : expicj) :, the cutoff x c can be elimi- 
nated by writing ( |B5[ ) using a normal ordered exponen- 
tial, tp satisfies (ip(x, r)ip(0, 0)) = [27t(«f t + The 
dynamics is governed by the Hamiltonian 



If = / d.v''l_ : (d x cj>) 2 



(B6) 



To make contact with the CFT notation, we first ob- 



serve that the normalization of ip in (B5) differs by a 



factor of V27T from (Bl|. Consider a finite system with 



periodic boundary conditions of length L, so that vpr+ix 
is defined on a cylinder, and introduce the radial variable 



z(x,t) 



a 27r(vFT-\-ix) j L 



(B7) 



The fermion field on the cylinder then has the form sim- 
ilar to piL 



2nz 



z— z(x,t) ■ 



(B8) 



Aside from the \/2ir difference in the normalization, this 



is equivalent to (B5 ) for L — >• oo. 



The Hamiltonian ( B6 ) corresponds to the lowest mode 



of the energy momentum tensor in the radial quantiza- 
tion. 



with 



Ln = 



2m 



dzzT(z) 



(B9) 



(BIO) 



where the integral is on a circle \z\ = e 27r " FT / L . It can 



readily be seen that (|B7j) , (|B9j) and (|B10|) imply that 
13 1 and (B6| are equivalent. 

Now consider the k fields <p a , along with (j) u — o a consid- 



ered in Section IV. A, which are related by (4.13 1. (Again, 
we consider only a single (ie right moving) chiral sector, 
and omit the superscript R). In the notation defined 
above, these satisfy 



(<t>a(z)<f>b(w)} 
(<Mz)<Mw)> 



Sab^og(z - w) 
-<W log(z - w 



and the energy momentum tensor is 

T [SU(2) 1 ]" = ~\Yj 
a 

l 

~ ~2 

For each of the k channels, the operators 

Ja = id z4>a 

J± = J m „ ± iJ v n = V2e ±l ^ 



(Bll) 

(B12) 
(B13) 

(B14) 



form an SU (2)i current algebra. Using the fact that 
: 3 a ■ 3 a : = —3 : {d z 4> a ) 2 :, we may write the energy 



momentum tensor in a way that reflects the \SU (2)i] k 
symmetry, 



T, 



[SU(2) 



(B15) 



The diagonal subalgebra generated by J = ^ a J„ 
forms a SU(2)k current algebra. The corresponding en- 
ergy momentum tensor will be proportional to J • J (we 
now suppress the normal ordering symbols, for brevity). 
The coefficient can be deduced by using (B14) to com- 
pute 



J J 



(k + 2)(d z 4> p ) 2 -2(d z ^f 



2 E e 

a=tb 



V2(d a 



(B16) 



If we require that the (d z 4> p ) 2 terms in (B13l and (B16) 
coincide, then an expression similar to (B4) shows that 
J z = 



2kd z 



bp has the appropriate scaling dimension 
A = 1. We then recover the Sugawara energy momentum 
tensor, 

1 



T, 



SU(2) k 



2(fc + 2) 



J J 



(B17) 



T. 



It follows that we may express T^ su ^ 2 ) 1 ] k m (B13) as 



SU(2 ) 



(B17), and T, 



T [SU{2) 1 ] k /su(2) k , where T S u(2) 



[SU(2) 1 ]"/SU(2) k 



is given in 
is the remainder. More- 



over, Tsf7(2) fe may be further decomposed into T ( 



T. 



SU(2) k /U(l), 



(7(1) -f 

where the U(l) term, generated by J z 



is simply the {d z 4> p ) 2 term in ( |B16[ ), and T S u{2) h /u(i) 



i s the rest. T his leads to the final results quoted in 
( |4.27|4.28|4.29 |. 



Appendix C: Operator Product Expansions 

In this appendix we provide the details of the calcula- 
tions for Eqs. 4.30| - [4~33| that show that the operators 
VP and T defined at the end of Section IV. A are pri- 
mary fields of the SU(2)k/U(l) Zk parafermion sector 
and the [SU(2)i] k /SU(2) k sectors respectively. The key 
is to compute the singular terms in the operator product 
expansion of the energy momentum tensors for the coset 



sectors (given in Eqs. 4.27 4.28 4.29 and discussed in 



the previous Appendix) with these operators. 

The necessary terms involve two kinds of products, 
which it is useful to consider separately. First, using 
Wick's theorem and (B2|, the OPE of the quadratic 



terms in T with an exponential operator takes the form 
- \{d z $ a {z)) 2 e* 5 -^ 

= { \D\ 2 /2 iD-dJ a (z) \ iD.$„( w ) ,py> 
\ (z — w) 2 z—w J \ / 



_ ( \D\ 2 /2 

(z—w) 2 



iD-4> a (w) _|_ 



For brevity we have suppressed the normal ordering sym- 
bols. This shows that the operator e 



iD-<f> a 



is a pri- 



mary field of the [SU{2) 1 ] k theory (described by (B13)) 
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with scaling dimension A = \D\ 2 /2. In particular, us- 
ing (4.11), this shows that A$ = \d a \ 2 = 1 — 1/fc and 
A T = A H = |d a - ( 4| 2 /4-l/2. 



Using the Baker Hausdorf formula and ( B2 ) , we may 
show that 



a iD 1 -tj> a (z) iD 2 -<i> a O) 



(z-w)- D i- B 2 



(C2) 



l+i(z-w)D 1 -d lll 4> tT {w) i(Di+D ? .)-$„(w) , 

(Z- W )- D 1 S 2 



The OPE's for * involve (fcTl) with D = V2d a , along 



with ((C2]) with D 1 = y/2(d a - d b ) and D 2 = V%d c . In 



this case, D\ ■ D 2 = 2(S ac — <5f, c ), so there are singular- 
terms in t he O PE only for c = b ^ a. Using the fact 
(from Eq. 4.9) that ^2 b ^ a d a — d b — kd a , we then find 



^2 e iV2(Za-z b )-$ ff (z) e 



i</2d c -<p a (w) 




Combining ([Cl]) (with \D\ 2 / 2 = 1 - 1/fc) and ( |C3| ), leads 
immediately to the results (4.30) and (4.31) quoted in 



Section IV. B. That Tym^ = is obvious because ^/ 
does not depend on <j) p . 

Thus, we have established that "J is a primary field 
of the SU(2)k/U(l) sector. Vl 1 are a bosonized repre- 
sentation for Z k parafermion operators^. They can be 



combined to define more general operators, 



*l = A k>l 



(C4) 



a 1 <...<a l 



For an appropriate normalization A 



k.i 



y/l\(k-l)\/k\ 

these can be shown using (|C2j) to satisfy the OPE's for Z k 
parafermions discovered by Zamolodchikov and Fateev^. 

The OPE's for T and S involve ([ClJ with D = (d a - 
d b )/V2 (\D\ 2 = 1), along with (|C2j> with D 1 = V2(d a - 
db) and D 2 = (d c — rfd)/v2. It follows that D\ ■ D 2 = 
<5ac + Sbd — 5 a d — 3bc- This leads to a l/(z — w) 2 singularity 
for a — d and b — c. In addition, there is a l/(z — w) 
singularity for a = d ^ b ^ c or b = c ^ a ^ d. After an 
analysis similar to ( C3 ) it follows that 



E < 

a^b,c^d 



-d b )-(fi^(z) e ±i(d 



2d u 



(z — w) 2 Z - 



E 

c^a^b 



2 cos 



((24 



a^b 

- d„. 



;)-0„/v^ C5 ) 
Z<r/V2) 



(z — w) 

Since the last term is independent of the sign in the expo- 
nent, it cancels for T. Combining (C5) w ith (|C1[ ) (with 
\D\ 2 /2 = 1/2 then leads to (|4~32]) and (|4.33[). Again, 



the U(l) term vanishes because T is independent of <fi p . 
Thus, T is a primary field of the [S , C/(2)i] fc /S'C/(2) fc sec- 
tor. 

For S, the last term does not cancel. The OPEs 



for both T. 



SU(2) k /U(t 



and Tr 



[SU(2) 1 ] k /SU(2) k 



are both non 



zero and do not have the form of a primary field. Pre- 
sumably, 5 can be written as a sum of terms that 
are products of primary fields in the SU(2)k/U(l) and 
[SU(2) 1 ] k /SU(2) k sectors. 
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